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ABSTRACT 

We use I band imaging to perform a variability survey of the 13Myr-old cluster 
h Per. We find a significant fraction of the cluster members to be variable. Most im- 
portantly, we find that variable members lie almost entirely on the convective side of 
the gap in the cluster sequence between fully convective stars and those which have 
a radiative core. This result is consistent with a scenario in which the magnetic field 
changes topology when the star changes from being fully convective, to one containing 
a radiative core. When the star is convective the magnetic field appears dominated by 
large-scale structures, resulting in global-size spots that drive the observed variability. 
For those stars with radiative cores we observe a marked absence of variability due to 
spots, which suggests a switch to a magnetic field dominated by smaller-scale struc- 
tures, resulting in many smaller spots and thus less apparent variability. This implies 
that wide field variability surveys may only be sensitive to fully convective stars. On 
the one hand this reduces the chances of picking out young groups (since the convec- 
tive stars are the lower mass and therefore fainter objects), but conversely the absolute 
magnitude of the head of the convective sequence provides a straightforward measure 
of age for those groups which are discovered. 
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1 INTRODUCTION 

Despite the widespread use of the photometric variabil- 
ity of young stars to study their rotational evolution (see 
iHerbst et al.ll2007l . for a review), its underlying causes re- 
main poorly understood. The variability itself appears to be 
due to a combination of (in order of i mportance) cool spots , 
hot spots, obscuration and pulsation fe xzesiriski et al I I 19991 : 
ICarpenter et al.ll200ll ; iGrankin etldlboOTl, bOOST). We know 
these stars have kG magnetic fields (e.g. |johns-Krulill2007l ). 
and all these sources of variability are probably related to 
those fields. The cool spots are large starspots, the hotspots 
are probably due to accretion funnelled onto a small region 
of the star by the magnetic field, and the obscuration is again 
due to magnetically confined accretion structures. A largely 
unexploited route to elucidate the underlying physics of the 
resulting optical variability is to study its characteristics as a 
function of fundamental stellar parameters such as magnetic 
field, age and accretion rate. Variability caused by pulsation 
will be of order millimagnitudes and so difficult to observe. 
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Additionally, variability caused by obscuration due to cir- 
cumstcllar, or disc, mater ial is expected to be negligible for 
h Pe r, as at » 13 Myrs (jSlesnick et alj|2002l ; IMavne et ail 
l2007t ) most of the associated discs and circumstellar ma- 
terial would probably have dissipated. However, variability 
caused by large cool spots or accretion hotspots generates 
significant magnitude changes over periods of hours to days, 
and therefore is relatively simple to detect. Here we study 
the degree of variability, caused by hot and cool spots, as 
a function of stellar mass in the young cluster h Per. The 
mass range covered by our data includes stars either side of 
the transition from low-mass fully convective to stars with 
radiative cor e s (the point in the colour-magnitude diagram 
IMavne et al.l (|2007f ) call the radiative-convective gap). 

The target of our investigation, h Per, and its com- 
panion x P er are y° un g) bright open clusters that have 
been the subject of detailed observations for s eventy years 
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.Crawford et all 1 197(1 iTapia et all 
Waelkens et alJll990l ). More recently, attention has 
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been focused on their mass functions and m ass segregation 
jSlesnick et al.ll2002l : iBragg fc Kenvonll2005l ). We adopt an 
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age of 13 Myr (|Slesnick et alj|2002l ; iMavne et al.ll2007l ), al- 
though the precise value is not crucial for our work. 

The remainder of this paper is laid out as follows. In 
Section[2]we describe the data acquisition and reduction pro- 
cesses, and the selection criteria used to identify variability. 
In Section[3]we present a colour-magnitude diagram (CMD) 
of the variable stars, alongside a brief statistical analysis of 
the gap detection significance. In Section [4] we discuss our 
findings and their implications. Finally, Section [5] is a short 
conclusion. 



2 METHOD 

The observations were taken using a Sloan i filter with the 
2.5m Isaac Newton Telescope (INT) on La Palma, using the 
Wide Field Camera with four EEV CCDs. The run spanned 
16 nights, from the 22nd September to the 7th of October, 
2004. h Per was observed on 12 of these nights, for around 
2.5 hours per night. Exposures of 5, 30 and 300 seconds were 
obtained. There were 213 good frames taken with 300s expo- 
sure, and 110 good frames for each of 5 s and 30 s exposure. 



2.1 Data reduction 

The data were reduced, and then photo metry extracte d us- 
ing the optimal ext raction algorithm of iNavlorl (|l998l ) and 
Navlor et al.l (l2002h. with the modifications introduced by 
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(2005). The images were bias subtracted us- 
ing a single master bias frame, corrected for flat-field effects 
using twilight sky fiat-field images, and defringed. For each 
exposure time we chose a single frame, and then searched 
the combination of the three to create a master catalogue. 
The objects were then matched to a 2MASS catalogue of 
our field to produce an astrometric solution, with a mean 
rms discrepency in position of 0.1 arcseconds. 

For each image positional and rotational offsets were 
determined with respect to the master catalogue. Using the 
list of candidate objects we performed optimal photometry 
on individual frames to produce fluxes for each star in each 
image. 

Next a "profile correction", the logical analogue of an 
aperture correction in standard aperture photometry, was 
applied to correct the flux measured for our objects to a 
somewhat larger aperture (although the crowded nature 
of the field limits this). A number of bright, unsaturated 
stars in each CCD provided the profile correction for that 
CCD, which is a polynomial function of position. Although 
our photometry is relative, and not tied to any particu- 
lar system, this correction is still important because the 
spatial variation of the PSF is a function of the seeing 
ILittlefair et al.ll2005T l. 

From this point we separated the photometry by ex- 
posure time into three datasets. For each dataset a rela- 
tive transparency correction was applied to all the result- 
ing photometric measurements, to normalise any differences 
between frames arising from variations in airmass or trans- 
parency. This was achieved by selecting a subset of relatively 
bright, unvarying stars and using their average magnitude to 
define a reference brightness for each frame. Th i s was done 
in the same way as described in iLittlefair et~ai] (|2005l ). 



The photometry for a star in each dataset is only rel- 
ative photometry with respect to the other stars in that 
dataset. To produce the final catalogue, the datasets were 
renormalised using the magnitude ranges for which they 
overlapped. We then applied an arbitrary zero point to our 
instrumental magnitude so that our observations would be 
in a natural magnitude system approximately comparable 
with Sloan i. We denote magnitudes in this system by the 
symbol ii . Of course many stars have magnitudes from more 
than one dataset, but for consistency it is important we cre- 
ate each lightcurve from a single dataset, and that we use 
the same dataset for stars of similar mean magnitudes. Thus 
the 5 s dataset was used for all stars brighter than an ii of 
15.2. Lightcurves for stars with 15.2 < ij < 17.5 were taken 
from the 30 s dataset, while for stars fainter than ii = 17.5 
lightcurves were taken from the 300 s dataset. These char- 
acteristic or limiting magnitudes are shown plotted over the 
data as dashed lines in Figure [T] 

2.2 Variable star selection 

The identification of variable stars involved a set of distinct 
selection criteria. First, the dataset was stripped of stars 
flagged as problematic during the data reduction process. 
Data quality checks made during this process included iden- 
tifying duplicate detections of the same star, stars contain- 
ing bad or saturated pixels, pixels with counts in the non- 
linear response range of the detector, and pixels with neg- 
ative counts. Additionally, stars possessing a PSF that was 
not point-like were flagged for non-stellarity. A number of 
spurious variables at the corner of CCD 3 were removed due 
to high vignetting, using a simple coordinate cut. We then 
selected only those lightcurves with a mean signal-to-noise 
per frame of more than 10. 

Having dealt with the most obviously problematic ob- 
jects, we next considered the lightcurves themselves. The 
ideal datset would have the same number of datapoints for 
each lightcurve, allowing us to use a simple \ 2 cut t o identify 
variable stars. Of course many lightcurves have missing dat- 
apoints for a variety of reasons, which means that we should 
choose a different value of \ 2 for each lightcurve if we require 
identical probabilities that the lightcurves show variability. 
In practice, for a large number of degrees of freedom the 
distribution of reduced % 2 is such that the probability as- 
sociated with a given xl is almost independent of v. This 
means if we pick only those lightcurves with a large num- 
ber of datapoints, we can apply a simple xt cutoff. Such 
an approach has the advantage that even if our model of 
the uncertainties is not perfect (in particular there may be 
correlated noise in the data), the cut-off will be consistent. 
So we required a minimum of 60 good datapoints from the 
5 and 30 second datasets (from a possible 110), and 120 
good datapoints from the 300 second dataset (from a pos- 
sible 213). To make our analysis more robust against sin- 
gle spurious datapoints we discarded the datapoint with the 
highest individual \ 2 with respect to the weighted mean, 
before re-calc ulating the mean , and hence \ 2 - 

Following ILittlefair et all (|2005T > . we chose a threshold 
of xl > 10 as a cutoff indicating significant variability in 
the 300 s dataset. The correct \ 2 limit for the 30 s dataset 
was then determined by considering stars falling within the 
overlap region between the two datasets. Varying the value 
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Figure 1. CMD of the variable population, after combining the 
5 s, 30s and 300s datasets (724 objects). The dashed lines indi- 
cate the limiting magnitudes (discussed in Section 12.11 used to 
combine the datasets. The fact that the sequence appears unbro- 
ken across the boundaries means that the selection criteria will 
not be a function of magnitude and therefore sensitivity. 



of xt applied to the 30 s dataset, we found that a threshold 
°f xt > 3 recovered a number of variable stars in this region 
that was comparable to the number of stars detected in the 
300 s dataset. Similarly, by considering variables identified 
in both the 5 s and 30 s dataset, it was found that a xt > 3 
threshold applied to the 5 s dataset recovered all variables to 
ii ~ 15.2. This process yielded 724 significantly varying ob- 
jects from the initial list of 19860 objects with good quality 
lightcurves. 



2.3 Periodic Variability 

One might hope to reduce the contamination still further 
by considering only those stars which have periodic modu- 
lation. Although our dataset is not well suited to identifying 
periods, caused by magnetically formed spots on the stellar 
surface, in the likely range for PMS stars (a few hours to 
«20 days), we did carry out an exploratory period search, 
summarised in Table [1] The perio d search was pe r forme d 
following the method described in iLittlefair et al.l {2005). 
Final analysis of the shortlist of variables was made by vi- 
sual inspection of individual lightcurves. Here we show our 
best candidates for periodic variability, along with a note of 
their position in the CMD. The number of periodic variables 
is clearly too small for our analysis, and so we instead use 
the catalogue of variable stars to create colour-magnitude 
diagrams. 



2.4 Colour-magnitude diagrams 

To create CMDs for the variable stars we cross matched our 
ca talogue of var i able s tars against the h and \ P er catalogue 
of iMavne etaH i|2007h . This gave the V vs V-I CMD with 
724 objects shown in Fig.[T] To enhance the contrast between 
the cluster and the contaminating stars (see Section [3]) we 
also cre ated CMDs f r om b oth our variable star catalogue 
and the IMavne et ail (|2007l ) catalogues for a 10 arcminute 




Figure 2. CMD of photometry towards h Per showing a clear 
sequence of probab le members, including variables, taken from 
IMavne et al.l J2007T) (upper panel), and variables only (lower 
panel) after a 10" radial cut centred on h Per. The upper panel 
has 3840 objects; the lower, 230. The number of variables drops 
at V-I < 1.5. 



radius region centered on the cluster (a = 02 18 58.76, S — 
+57 08 16.5 J2000). With this spatial cut, these catalogues 
have 230 and 3840 objects, respectively, plotted in Fig. [2] 



3 RESULTS 

Fig. [T] shows a clear pre- main-sequence population running 
diagonally from approximately 14th to 20th magnitude in 
V . However, background contamination is also evident. The 
contamination is more obvious in the upper panel of Fig. [2] 
where it covers a large, wedge-shaped area that cuts laterally 
across the sequence. The majority of this wedge is likely to 
be dwarfs, but the additional "finger" at V-I ~ 1.6 is prob- 
ably due to background giants. As described in Section 12.41 
to reduce this contamination we selected only those vari- 
able stars within 10 arcmin of the cluster core. Fig. [2] (lower 
panel) shows the variables remaining after this process. We 
can clearly see that the density of stars falls dramatically as 
we progress up the sequence. This is quite different to the 
distribution of cluster members (Fig. [2] top panel), which 
shows a narrow sequence up to around 15th magnitude. 



3.1 Histogram analysis 

Although it is clear from visual inspection of the sequence 
that the density of variable stars in the CMD falls dramati- 
cally at brighter magnitudes, additional analysis is required 
to investigate this behaviour in a robust way. Because we 
wish to examine the distribution of variables along the clus- 
ter sequence, we construct histograms according to the fol- 
lowing scheme. A straight line was fitted to the central por- 
tion of the sequence, running from 1.0 to 2.0 in V-I. Two 
lines perpendicular to this sequence line in V-V—I space (i.e. 
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ID 


RA 


Dec 


V 


V uncertainty 


V-I 


V-I uncertainty 


Period (days) 


1.01 1293 


02 20 44.582 


+56 54 56.35 


18.948 


0.013 


1.935 


0.017 


0.1186993 


1.02 1325 


02 17 46.554 


+57 04 18.08 


18.312 


0.012 


1.337 


0.019 


0.1464187 


1.02 2157 


02 18 01.071 


+56 53 56.73 


18.770 


0.015 


1.159 


0.024 


0.2099175 


1.04 0251 


02 20 21.400 


+57 03 22.03 


18.685 


0.010 


1.994 


0.014 


0.1160178 


1.04 1906 


02 20 28.387 


+57 03 03.66 


19.325 


0.015 


1.854 


0.021 


0.3029911 



Table 1. Periods for the five strongest periodic variable candidates in our sample. 




Bright cut 



Figure 3. Placement of cuts along the line of the sequence of 
h Per variables. Stars falling within either of the cuts are repre- 
sented as crosses. Identical cuts were used to analyse the general 
cluster membership. 



appearing perpendicular to the sequence when both axes are 
plotted with units of equal size) were then taken, a bright 
cut crossing the sequence at 1.0 in V-I, and a faint cut 
crossing at 1.66 (see Figure [3}- We selected the stars within 
0.4 mags of these lines and then pulled them back onto the 
cut centre-lines and binned the m in V-I. Iden t ical cu ts were 
used for the h Per catalogue of iMavne et al. | \200H ). These 
cuts produce the histograms of Figure [3] 

The histograms show in detail what can be seen qual- 
itatively by eye in the CMDs of Fig. [2] While the general 
cluster membership continues to relatively blue colours, the 
variability dies out blueward of about V-I = 1.5. In the 
bottom-left panel of Fig. [4] which shows the faint magni- 
tude (lower) cut, the variables clearly trace the sequence at 
1.6—1.7 in colour. The sequence can also be seen in the his- 
togram for the general cluster (Fig. [4] top-left), as a shoul- 
der of excess sources over the background, although in this 
case the contamination dominates between 1.0-1.5 in V-I. 
That these stars are mostly contaminants is evidenced by 
the large number of background objects present in the sky 
around the cluster at these colours (dashed line, normalised 
to the same area as the selected area of the cluster). Note 
that the blue end of the contamination appears to be shifted 
bluewards, relative to the for eground cluster. This coul d be 
due to differential reddening (|Capilla fc Fabregatll2002l ). 

In the brighter magnitude cut (Fig.[4] top-right) the 
cluster is significantly less contaminated, with more than 
half of the stars additional to the background. However, 
looking at the variables (Fig. [4] bottom-right) presents a 
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Figure 4. Histograms of the lower cut (left hand figures) and up- 
per cut (right hand figures), binned along the line of the cut. The 
upper panels show non- variables, while the lower panels present 
only stars exhibiting significant variability. The solid line indicates 
objects falling within a 10' circle of the h Per centre. The dashed 
line is a count of the stars falling outside that cut, normalised 
to the same area as the cut. The dashed lines thus measure the 
background contamination present in the field. 



different picture. There are relatively few variables present 
either in the background or in the sequence. In fact, the 
number of variables is consistent with a scenario in which 
there are no variables present in the sequence. 

The position of the two cuts was chosen to bracket the 
ra diative-convective g ap (as defined and identified in h Per 
bv lMavne et al.ll2007l ). where the stars move from the fully 
convective regime to development of a radiative core. There- 
fore, these data and the chosen cuts allow us to analyse 
changes in the fraction of variables either side of this gap. In- 
tegrating under the histograms of Figure [4] allows us to con- 
struct the background-subtracted numbers of variable and 
non-variable stars, and calculate the fraction of stars which 
are variable, both above and below the gap. We find that 
below the gap 37.9 percent of stars are variables, above it 
—1.7 percent. The latter is negative as the total (background 
plus cluster) number of variables is so small that its associ- 
ated noise means that background subtraction can push it 
below zero. This means that we actually see less variables 
in this region of the CMD than the mean number we would 
expect to be contributed by background objects. To test the 
significance of this difference we formed the null hypothesis 
that the fraction of variables above the gap is the same as 
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that below the gap, i.e. 37.9 percent. We then used the Pois- 
son distribution to simulate observations of the stars above 
the gap, using parent populations identical to our observa- 
tions, save that for cluster stars above the gap we changed 
the fraction of variables to 37.9 percent. In one million sim- 
ulations we never obtained an observed fraction of variable 
stars less than or equal to -1.7 percent. This allows us to 
reject the null hypothesis that the ratio of members to vari- 
able stars is constant either side of the RC gap with great 
confidence, and conclude that the variable stars do indeed 
die out as one moves above the gap. 



4 DISCUSSION 

4.1 A break at the RC gap? 

IStolte et ail l|2004 ) and iMavne et all < |2007f > show that in 
CMDs there is a discontinuity in the stellar sequence be- 
tween pre-main-seq uence and main-sequence stars, which 
IMavne et al.l l|2007l ) term the radiative-convective (RC) gap. 
Below this gap theory predicts that the stars are fully con- 
vective, whilst above it they possess radiative cores. This 
gap closes with age, becoming very small for objects as old 
as h Per. For stars above the gap, one model postulates 
that the magnetic field is generated at the tachocline be- 
tween the radiative core and the convective envelope. In 
this model, it is believed that shearing processes in this layer 
generate a magnetic field, although the full pictur e may well 
be more complicated fe.g. IChoudhuri et alj|l99a ). For fully 
convective stars the source of the field is more controver- 
sial. On the one hand a turbu lent velocity field m ight create 
a sm all-scale magnetic field tourney et al,l fl993). Alterna- 
tively [ChabrieL&iKtikei feom showed that a stratified ro- 
tating medium might create a large-scale magnetic field (the 
a 2 dynamo). 



IMavne et al.l (|2007T ) identify the position of the RC gap 
in h Per as V = 17.66 and V-I = 1.33. It is clear from Fig.0 
that there is a marked increase in the number of variable 
stars below this position in the sequence. This supports the 
conclusion that the mechanisms giving rise to the variability 
are a function of convective processes in the star, and are 
curtailed by the development of a radiative core. 

One interpretation is that stars below the gap (which 
are fully convective) have large-scale spot structures, as evi- 
denced by their variability, implying equivalently large-scale 
magnetic fields. In contrast, stars above the gap, which have 
a radiative core and greatly reduced convective zone appear 
to be less dominated by large-scale magnetic field structures, 
as evidenced by the dramatic absence of variability arising 
from spot structures. This implies either an absence or suffi- 
cient uniformity of distribution of spots, resulting in a mod- 
ulation, if any, below our threshold of detection. Thus our 
observations are consistent with the idea that whilst fully 
convective stars have large-scale structure in their magnetic 
fields, once the stars develop a radiative core the spots be- 
come more evenly distributed over the surface. Whilst it 
can still be argued whether the magnetic fie ld is due to 
convective turnover l|Chabrier fe Kiikerl |2006[). turbulence 
|Durnev et alJll993h or fossil field |Tavlerlll987h . what this 
study makes clear is that the fully convective structure is 
a crucial ingredient for the maintenance of the large-scale 
magnetic fields of T Tauri stars. 



4.2 Comparison with older field stars 

Further support for a change in magnetic field topology at 
the RC gap can be found from recent studies using Zeeman 
Doppler imaging of older field dwarfs, where the same physi- 
cal mechanisms are apparent, i.e. fully convective lower mass 
stars and higher mass ojects with radiative cores. These 
studies have found significant evidence for a sharp reduction 
in magnetic flux, contained within large scale structures, as 
one moves across the boundary (in mass) betw een fully con- 
vectiv e sta rs and those with radiative cores. iDonati et al.l 
l|200Sl) and lMorin et al l (|2008T ) compare observations of the 
Stokes V line profiles for old field dwarfs about the boundary 
of radiative core pr esence (at M* ~ 0.3 5Afp) using the evo- 
lutionary models of IChabrier fc Ba raffc 199?]), to simulated 
profiles from estimated magnetic maps of the stellar sur- 
face. The Stokes V measurements are sensitive to large scale 
magnetic fields, rather than small scale structures, which 
hav e opposite polaritie s th at cancel in the St okes V pro- 
file. iDanatTeFall <|2008f > and lMorin et all (|2008l ) find a clear 
difference in magnetic fields, with large scale poloidal, ax- 
isymmetric fields dominant for fully convective stars, and 
only toriodal and non-axisymmetric poloidal fields apparent 
over large scales for stars with radiative cores. Measurements 
of X-ray activity and the total magnetic field from Zeeman 
broadening show little change at the rad iative core bound- 
ary (see discussion in IDonati et al.l 120081 ) , implying only a 
structural change in the magnetic field, possibly due to the 
theorized change in dynamo mechanism. Additionally, for 
fully convective stars these large scale structures were ob- 
served to be stable over « lyr with stability over only a few 
months for the magnetic structures of stars with radiative 
cores. 

Further work using Stokes I ob servations, which issen- 
sitive to the total magnetic field, inlReiners fc Basril (2009) 
has similarly observed a sharp change in apparent field struc- 
ture over the transiti on from fully convective stars to those 
with a radiative core. iReiners fc Basril (|2009h find, by com- 
paring Stokes V and I components, that across this bound- 
ary (i.e. the RC gap) the percentage of total flux stored in 
large-scale structures drops from 15 to 6%, with the presence 
of a radiative core. The mean total flux does no t show a com- 
parable jump. Interestingly, IReiners fc Basril (2009) show 
that 85 percent of the flux is not detected in Stokes V and 
is therefore is present in small scale strucures. Additionally, 
they note that the weak-field assu mption and least-sq uares 
deconvolution met hod applied in IDonati et al.l (|2008f ) and 
iMorin et al.l (|2008l ) could lead to large scale fields of sev- 
eral kG being missed. This may mean that the change in 
magnetic field structure with the transition from full convec- 
tive to radiative core could simply be due to an increase in 
field strength of the large-scale structures over the detectable 
threshold. This scenario is currently untestable, but is un- 
likely given our observations. Our data indicate a change 
in photometric variability and therefore magnetic field stru- 
cure with the presence of a radiative core. If this is the same 
mechanism as is postulated for the older field stars, and is 
indeed the result of a transition in stellar structure, then 
the change in variability observed across the RC gap arises 
from a fundamental change in magnetic field topology, itself 
a consequence of a change in the dynamo mechanism. 



6 Saunders et al. 



4.3 Implications for future variability surveys 

Variability will be an important technique for identify- 
ing young stellar associations in large area variability sur- 
veys such as the Large Synoptic Survey Telescope (e.g. 
iBricefio et al.ll2007f ). Fig.[2]shows what such a survey should 
discover. Only the sequence below the RC gap is visible in 
the variability data, which reduces the chances of finding 
the cluster. However the variability data clearly show where 
the RC gap begins. Since the position of the RC gap fades 
with age, measuring this position provides a new technique 
for measuring the age of a group, which will be useful well 
beyond 20 Myr. T his contrasts sharply with the results of 
iMavne et alj ( 20071 ) , who concluded that the distance across 
the RC gap would become unusable as an age indicator for 
groups older than h Per, since the gap closes to be extremely 
small. Thus it appears that although the RC gap in CMDs 
alone may no longer provide ages for these older systems, a 
CMD of the variable stars could do so. 



5 CONCLUSION 

In conclusion we have measured, and proved the statisti- 
cal significance of, a reduction in the number of variable 
stars, compared to non- variables, as one moves from below 
to above the RC gap in h Per. We have discussed the connec- 
tion between this reduction in variability acoss the gap and 
underlying changes in the stellar magnetic fields. We suggest 
that this distinct reduction in variability could be due to a 
change in the magnetic field topology caused by a change 
in stellar structure, from fully convective stars to those with 
radiative cores. This idea implies that the fully convective 
stars retain large-scale stable surface magnetic fields, which 
lead to long-lived cool spots on the stellar surface, and there- 
fore long-term variability. In contrast, for stars with a radia- 
tive core the magentic field becomes dominated by smaller 
scale structures producing more transient and less significant 
variability. It is clear that variability surveys will therefore 
be biased towards detection of the fainter, fully convective, 
stars within clusters. However, the bright edge of the RC 
gap, where the variability dies out, could be used as an age 
indicator. 
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